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We study spin-half and spin-one Heisenberg models in the limit where one dimensional (1-D) 
linear chains, with exchange constant Ji, are weakly coupled in an anisotropic triangular lattice 
geometry. Results are obtained by means of linked-cluster series expansions at zero temperature 
around different magnetically ordered phases. We study the non-colinear spiral phases that arise 
classically in the model and the colinear antiferromagnet that has been recently proposed for the 
spin-half model by Starykh and Balents using a Renormalization Group approach. We find that such 
phases can be stabilized in the spin-half model for arbitrarily small coupling between the chains. For 
vanishing coupling between the chains the energy of each phase must approach that of decoupled 
linear chains. With increasing inter-chain coupling, the non-colinear phase appears to have a lower 
energy in our calculations. For the spin-one chain, we find that there is a critical interchain coupling 
needed to overcome the Haldane gap. When spin-one chains are coupled in an unfrustrated manner, 
the critical coupling is very small (« O.OlJi) and agrees well with previous chain mean-field studies. 
When they are coupled in the frustrated triangular-lattice geometry, the critical coupling required 
to develop magnetic order is substantially larger (> 0.3Ji). The colinear phase is not obtained for 
the spin-one Heisenberg model. 



PACS numbers: Valid PACS appear here 
I. INTRODUCTION 

There has been considerable recent interest in the 
properties of two-dimensional (2-D) antiferromagnetic 
Heisenberg models. In the absence of frustration, the 
ground state phases and properties of these models are 
quite well understoodii^i^. On the contrary, a complete 
knowledge of the ground state phase diagram of frus- 
trated Heisenberg models is still lacking. 
In the present paper, wc further study the antiferro- 
magnetic Heisenberg model on the anisotropic triangular 
lattice^. For this class of models, the Hamiltonian can 
also be defined on a square lattice with nearest neigh- 
bor interaction J2 > and a second-neighbor interac- 
tion Ji = 1 along one of the diagonals of the squares, 
as shown in Fig. [TJ In the limit J2 — > 0, the model re- 
duces to 1-D decoupled spin chains along the diagonals. 
In the limit J2 = 1 the model becomes the Heisenberg 
model on the 2-D isotropic triangular lattice, for which 
there is strong numerical evidence of long-range antifer- 
romagnetic order— i^iii^. We are interested in studying the 
region < J2 < 1 to see how the order develops as the 
interchain coupling is increased. 

One reason for strong interest in these models is 
that the Heisenberg model with J2 ^ 1/3 provides 
the dominant terms in the Hamiltonian for the mate- 
rial Cs2CuCl4^. The exchange parameters of this mate- 
rial have been determined from high-field studie a^°'^^ , by 
measuring the excitation energies around the fully polar- 
ized limit. In zero field, these materials are found to have 
spiral long-range order in the ground state^^. The mate- 
rials also have non-zero Dzyloshinski-Moria (DM) inter- 
actions. Spin- wave theorie a^^i^^ can account for the spin- 
wave dispersion in the materials only after substantially 




FIG. 1: (Color online) Square lattice with coupling constant 
J2 (dashed black line) along the horizontal and vertical axis 
and Ji (solid red line) along the diagonal 



modifying the exchange parameters. On the other hand, 
series expansions around the spiral state give an excel- 
lent account of the main peaks observed experimentally 
in the spectra throughout the Brillouin zone, with the ob- 
served exchange parameter o^^'^^i^^ . In another approach 
to the spectra of these materials, Kohno et &\t^ have 
shown that many features of the experimental spectra, 
including the observed continuum can be well explained 
by considering weakly coupled Heisenberg chains, even 
though the interchain coupling is not too small in the 
material. 

In another recent study, directly relevant to the present 
work, Starykh and Balcntsi^ have considered the frus- 
trated weakly coupled chain problem using a Renor- 
malization Group approach. Their striking prediction 
is that a colinear-phase is stabilized in place of the 
classical spiral phase for sufficiently weak coupling be- 
tween the chains. The phase corresponds precisely to 
the four-sublattice phase that is known to occur in the 
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square-lattice Ji — J2 model at large J-^^iSL. Other 
analytical, numerical and variational studies have also 
been used to study this anisotropic triangular-lattice 
model^^i^^i^iiS^i^^. Several of these have found evidence 
for disordered spin-liquid phases in the model at weak 
coupling between the chains. Most notably the vari- 
ational calculations of Sorella and coworkers find two 
distinct spin-liquid phases in the model. Furthermore, 
the DMRG studies^l of Wcng et al. showed a very 
rapid exponential decrease in correlations perpendicular 
to the chains even for rather large interchain couplings. 
One possible concern with these studies is that periodic 
boundary conditions were used and that can play a role in 
destabilizing non-colinear phases if they occur at incom- 
mensurate wavevectors. Earlier series expansion studies^ 
also found that the energies from spiral-phase series ex- 
pansions and dimer expansions were very close. In gen- 
eral, numerical studies of weakly coupled chains in frus- 
trated geometries have been a challenging proble m^^i^^ . 

In this paper, we revisit this model for the spin-half 
case and also study the corresponding spin-one model. 
To our knowledge, this is the first quantitative numeri- 
cal study of the spin-one model. We use an Ising-type 
linked-cluster expansion method^ at zero temperature 
around different phases. Ground state energy and sublat- 
tice magnetization have been calculated for each phase. 
The knowledge of the exact behavior of the 1-D spin-half 
modcl^i, with its essential singularity in the energy and 
magnetization functions^, is used to improve the series 
extrapolations in the 1-D limit, and they are also used 
to get more accurate estimates of the ground state en- 
ergies with inter-chain couplings. While both colinear 
and non-colinear phases can be stabilized in our studies 
for the spin-half model, we always find the energy of the 
non-colinear phase to be lower. However, the results are 
quite sensitive to the way the scries are analyzed, espe- 
cially in the limit of weak interachain couplings, and this 
implies some uncertainties in our results that cannot be 
addressed by series expansions alone. Further numerical 
studies of these models would be useful. 

We have also investigated the spin-one model, with 
unfrustrated (square-lattice geometry) and frustrated 
(triangular-lattice geometry) interchain couplings. In 
this case, the 1-D limit corresponds to a Haldane gap 
phase. The Ising expansions are known to break down 
before the Heisenberg symmetry is restored, with a criti- 
cal point which is in the universality class of the 2-D Ising 
model. We find that when the chains are coupled in an 
unfrustrated manner a rather small interchain coupling 
{J2/J1 < 0.01) leads to Nccl order. On the other hand, 
in the frustrated geometry a much larger interchain cou- 
pling (J2/J1 > 0.3) is needed to obtain long-range order. 
In the latter case, we only find the spiral phase to be 
stabilized for Heisenberg models. The colinear-phase be- 
comes less and less stable with increased inter-chain cou- 
pling. It should be stressed that the analysis of Starykh 
and Balentsi^ was special for the spin-half case and hence 
there is no apriori reason to expect a colinear phase in 



the spin-one model. 

Our results are organized as follows. In Section [III we 
discuss the methods of series expansions. In Section IIIII 
ground state energy and sublattice magnetization for the 
spin-half chain are presented. In Section IIVI we discuss 
our results for the spin-half model on the anisotropic tri- 
angular lattice. In Section |V] the study of the spin one 
model is presented. Finally, in Section IVll we present our 
conclusions. 



II. SERIES EXPANSION 

The antifcrromagnctic Heisenberg model is defined by 
the Hamiltonian 

i/ = Ji5^S..S, + J2 5] S,.S,. (1) 

[hj] <»J> 

Here, refers to one of the diagonal next-nearest- 

neighbor pairs on the square-lattice shown in Fig. 1, with 
corresponding coupling constant Ji, while < i,j > are 
pairs of nearest-neighbors with coupling constant J2 . We 
set Ji = 1, and vary J2 in the range < J2 < 1- In 
the limit J2 = the model is equivalent to decoupled 
antiferromagnetic spin chains and is exactly solvable for 
spin-halSi. For J2 = 1, the model is equivalent to the 
Heisenberg model on an isotropic triangular lattice. In 
this limit the classical ground state has a 3-sublattice 
'120-dcgrce' order, which can also be regarded as a non 
colinear spiral with wave wector q — cos~^(— i) = 
The predicted CAF phase of Starykh and Balcntsi^ for 
small J2 is shown in Fig. [2^. In this phase, the spins are 
aligned antiferromagnetically along the diagonals and the 
vertical axis of the square lattice, and fcrromagnetically 
along the horizontal axis. A sketch of the classical spiral 
phase is shown in Fig. [2]d. It was found in earlier series 
expansion studies^ that away from the triangular lattice 
limit, quantum fiuctuations renormalizc the angle q with 
respect to the classical value. We will refer to this renor- 
malized spiral phase as the non colinear antiferromagntic 
phase (NCAF). 

To obtain a T = expansion about the CAF phase, 
the Hamiltonian is written as 

H = Ho + X{Hi + H2) (2) 

where 

Ho = JiY, S-S] + ^tS- (3a) 

<*.i> 

H,=J,Y.{StS-+S-S+) (3b) 

H2 = J2 J2 ^^tS~ + S-S+) (3c) 

<i,j> 
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(a) 




(b) 



FIG. 2: (Color online) (a) CAF phase: spins are aligned an- 
tiferromagnetically along the diagonals and the vertical axis 
of the square and ferromagnetically along the horizontal axis, 
(b) Classical spiral phase. 



and A is the expansion parameter. The limits A = and 
A = 1 correspond to the Ising model and the isotropic 
Heisenberg model, respectively. Hq is taken as the un- 
perturbed Hamiltonian while H' ~ Hi + H2 is the per- 
turbation operator. To obtain a T = expansion about 
the NCAF phase, we rotate all the spins so as to have a 
ferromagnetic ground state. In this case the Hamiltonian 
form in (l2]) is still valid but now 



Ho = Ji cos(2(7) S^S^ + 

+ J2C0s{q) S^S^ 



(4a) 



<hj> 



Hi=JiJ2 ^f^^ + cos(2(7)S'f S'^ + 



(4b) 



<hj> 



-sin(g)(5f5f-Sr5n 



(4c) 



where q is the wave vector of the NCAF phase. Once 
the ground state phase has been chosen, perturbation 
theory can be applied, leading to an expansion of H up 
to desired order in A for the ground state wave function 
(IV'GS >); the ground state energy and other correlation 



functions. The sublattice magnetization is given by. 



< M >= 



(5) 



The details behind the technique are discussed elsewhere 
and will not be repeated here. For a complete review see 
Ref^i^ 



III. ONE DIMENSIONAL SPIN-HALF 
HEISENBERG MODEL (J2 = 0) 

The ground state properties of the 1-D spin-half 
Heisenberg model at T = are well known from the ex- 
act solutions^. It is also known^ that the ground state 
energy and sublattice magnetization, as a function of A, 
have essential singularities of the form 



exp(- 



1 



(6) 



We first study the 1-D case, to see how well the series 
expansion methods can reproduce the exact results. Se- 
ries coefficients for the sublattice magnetization of this 
model are generated up to order 10 in A. Their anaylsis 
is carried out in two different ways: (i) using Integrated 
Differential Approximants (IDA) on the series obtained; 
(ii) using Biassed Integrated Differential Approximants 
(BIDA) on the natural logarithm of the same series. The 
first approach is the most straightforward and simply fits 
the known coefficients of the series to a homogeneous or 
inhomogeneous differential equation of the form 

^^'^^^S + ^^^^^^ + ^''^^^-^ + ^^^^^ " ° 

where Pk,Ql, Rm, St are polynomials of degree 
K, L, M, T respectively. The results reported in this sec- 
tion are obtained by setting the polynomial Pk and St 
to zero, which is equivalent to a Dlog Fade analysis 



ldf_ 
fdX 



RmW 
' Ql{\) 



(8) 



The second analysis method proceeds by taking the nat- 
ural logarithm of the calculated series in order to reduce 
the essential singularity © into an algebraic one 



ln(e 



1 



(9) 



with critical exponents 7 = — 0.5. A sigularity of this 
form can easily be analyzed by using an IDA type of 
analysis. Moreover, a better convergence can be achieved 
by biassing the exponent in the analysis. This simply 
means that the function / in ([8]) is forced to have an al- 
gebraic singularity of the form ^ with critical exponent 
7 = —0.5, as known from exact result. Fig. [3] shows the 
comparison between the two analysis methods and the 
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FIG. 3: (Color online) (a) Sublattice magnetization, normal- 
ized to unity in the Ising limit. The red crosses represent the 
approximants obtained by analyzing the natural log of the 
series with BID A. The grey triangles represent the approxi- 
mants obtained by analyzing the original series by standard 
IDA 



exact result. While the two methods give very similar 
results for A < 0.7, it is clear that they tend to disagree 
for A — > 1. The approximants obtained by Biassed IDA 
on the natural logarithm of the series, reproduces the 
behavior of the series close to the critical point much 
better. This is expected as it is extremely difficult to 
account for an essential singularity of the form ^ with 
a short, finite series without the biassing. On the other 
hand. Fig. [3] shows that, once the essential singularity 
has been reduced to an algebraic one, by taking the nat- 
ural logarithm of the series, even a finite series up to ten 
terms can reproduce quite well the exact behavior of the 
system all the way to the critical point. 



IV. SPIN-HALF HEISENBERG MODEL ON 
ANISOTROPIC TRIANGULAR LATTICE 

The ground state energy for the Hamiltonian (1) has 
been computed up to order 10 in A for both NCAF and 
CAF phases for different values of the interchain cou- 
plings J2. For the NCAF phase, we consider a range of 
g-values and minimize the energy with respect to q. For 
each value of J2, we then calculate the ratio series 



i?j2=o(A) 



(10) 



where Ej^ (A) is the energy series calculated at a specific 
value of J2 and Ej^^q(X) is the energy scries computed 
for the 1-D model. The idea behind this is that, if the 
series for J2 7^ has apparent singularities as a function 
of A, as a consequence of being close to the 1-D limit, 
we can eliminate its effect by taking the ratio. This al- 
lows us to evaluate the energy with increased accuracy. 
Moreover, to improve convergence, an additional term is 



FIG. 4: (Color online) Energy for CAF (red cross) and NCAF 
(blue stars) phases. The inset shows a zoom in of the region 
around J2 = 0.1. 



added to the Hamiltonian (2) as in previous studies^i^i. 



(11) 



For A ^ 1, the limit we are interested in, this has no effect 
on the Hamiltonian. The amplitude of the convergence 
term t is generally taken as 1 . The extrapolated series ra- 
tio PU)) is multiplied by the exact result for the 1-D case, 
E = — ln(2) + i, to obtain the energies at different J2. 
Fig. m shows the obtained results. In this plot, the error 
bars reflect the spread of well-behaved approximants. At 
J2 = 1 the energy for the NCAF phase is centered around 
—0.5508, in good agreement with earlier studiesi^. For 
< J2 < 1, Fig. 3] shows that, in our calculations, the 
NCAF phase has a lower energy than the CAF phase for 
all values of J2. The inset shows that even for J2 = 0.1, 
the lowest data point taken in our calculation, the NCAF 
phase appears to have lower energy than the CAF phase. 
However, while this result is suggestive that the NCAF 
phase is the correct phase of the model, we should em- 
phasize that our results are sensitive to how the series 
are analyzed, especially near the decoupled chain limit. 
Hence, further numerical studies of this model, looking 
in an unbiassed way at short distance spin correlations, 
would be useful. 



A. Sublattice magnetization 

In this section, we study the sublattice magnetization 
for J2 7^ 0. Because the sublattice magnetization van- 
ishes for the 1-D model, it is not useful to consider the 
ratio of the sublattice magnetization series with that at 
J2 = 0. 

We have calculated series expansion coefficients for the 
sublattice magnetization of the NCAF and CAF phases 
to order 10 in A. Zheng et al4J^ had earlier calculated 
the sublattice magnetization for the classical spiral phase 
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for 0.25 < J2 ^ 1- These are shown in Fig. [5l In the tri- 
angular lattice limit, the magnetization is 0.19(2), and 
as we move towards the 1-D limit, it begins decreasing 
almost linearly to zero for J2 < 0.5. We have tried var- 
ious approaches to analyzing the magnetization series. 
They lead to similar results but with no improved con- 
vergence. The series analysis is even less reliable in the 
CAF phase. It is possible that the results at small J2 arc 
strongly influenced by the nearby essential singularity at 

J2=0. 

Motivated by the idea that the nearby singularity at 
J2 = may be influencing the series analysis, we take 
the following approach: We subtract off for each value of 
J2 7^ the series for the 1-D limit (J2 = 0). This should 
eliminate the influence of the essential singularities. We 
then analyze the series with standard IDA. The results 
obtained in the region 0.1 < J2 < 0.35 following this pro- 
cedure are shown in Fig. O (red triangles). This analysis 
leads to very small values of the magnetization, consis- 
tent with the exponentially small values expected from 
the work of Bocquet et al.— . Unfortunately this analysis 
does not smoothly connect with the results at large J2, 
so it is not clear how far in J2 it should be continued. 
If this scenario is correct, there maybe a sharp quantita- 
tive change of behavior between small and large J2 . We 
note that other groups have even suggested various phase 
transitions as a function of J^r^i^l. A similar analysis for 
the sublattice magnetization series for the CAF phase 
gives only negative values. 
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FIG. 5: (Color online) Sublattice magnetization for NCAF 
phase, normalized to a classical value of 0.5 calculated by 
Zheng et al.— (blue squares and gray triangle) and in the 
present work after subtracting the ID series (red triangles). 
See text for details, (inset) Zoom in of the small interchain 
coupling region showing the data points calculated in this 
work. 
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triangular-lattice. Ground state energy and sublattice 
magnetization have been calculated for CAF and NCAF 
phases as for the spin-half model. We begin by showing 
in Fig. [S]the wave vector q, in units of tt, for the non co- 
linear antiferromagnetic phase (NCAF) as a function of 
J2. For each value of J2, the energy has been minimized 
with respect to q and the value of q for which E = Emin 
has been plotted. The analysis has been carried out for 
the spin-i and the spin-1 Heisenberg model. The clas- 
sical result q = cos^^(— ^j-) is shown by a solid black 
line. The NCAF phase for the spin-1 system is closer 
to the classical solution than the spin-i system, as ex- 
pected. Nevertheless, evidence for quantum effects in the 
properties of the ground state are clearly evident. 




In this section we present calculations for the ground 
state properties of the spin-one model on the anisotropic 
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FIG. 6: (Color online) Comparison between the NCAF phase 
wave vector q (in units of tt) for the spin-1 and spin-^ 2-D 
Heisenberg model on anisotropic triangular lattice as a func- 
tion of J2. The classical result is also shown by a black solid 
line. 



A. Phase Diagram for anisotropic square and 
triangular lattice spin-one models 

In the 1-D limit, it is well known that the Heisenberg 
system is in the Haldane gap phase^^. Furthermore, Ising 
expansions lead to a critical point before the Heisenberg 
symmetry is restored (Ac < 1)^. Once, the couplings 
J2 are turned on, we would like to follow Ac as a func- 
tion of J2 to see when it reaches unity. This will tell us 
the critical J2 required to close the Haldane gap for the 
Heisenberg system. 

In Fig. [71 we show the value of the anisotropy param- 
eter Ac as a function of the couplings J2 obtained by a 
Dlog Fade analysis of the sublattice magnetization se- 
ries. One can see a clear difference between the NCAF 
and CAF phases. The CAF phase becomes less and less 
stable with interchain coupling and is never realized for 
the Heisenberg models (A equal to unity). In contrast, 
the spiral phase become more stable with interchain cou- 
pling and for some J2 the critical A reaches unity. It 
turns out that this way of studying the critical J2 is less 
accurate than an alternative method discussed below. 
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FIG. 7: (Color online) Ac vs. interchain coupling J2 for NCAF 
and CAF phase of the spin-1 Heisenberg model on the trian- 
gular lattice. 



in which 
estimates for J|''** 



and 7 are free parameters. This leads to 
of approximately 0.01 and 0.33 respec- 
tively with exponent values which are around a third, as 
expected for 3-D models^!. The value of J|"* « 0.01 
for the former is consistent with chain mean-field the- 
ory estimates^. It is clear that this kind of frustra- 
tion significantly extends the Haldane gap phase. Note 
that this is very different from a material like CsNiCla 
where a given chain has six neighboring chains arranged 
in a triangular-lattice geometry. In that case, the chains 
themselves remain unfrustrated and frustration only af- 
fects the relative spin-orientation between the chains^. 
The latter is much less effective in preventing long-range 
order. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
h 

FIG. 8: (Color online) (blue stars) Magnetization vs. inter- 
chain coupling J2 for the NCAF phase of the spin-1 Heisen- 
berg model on the triangular lattice, (grey triangles) Magne- 
tization vs. interchain coupling J2 for the Neel phase of the 
spin-1 Heisenberg model on the square lattice, (dotted lines) 
Fit of the magnetization curve to a power law of the form 
given in equation (|12|) 



A better way to compute critical J2 needed to close 
the Haldane gap is to consider the sublatticc magnetiza- 
tion as a function of J2 coming from the triangular-lattice 
side (large J2). The series analysis should be most accu- 
rate when the system is well ordered. In Fig. [51 we show 
the calculated sublattice magnetization for the Heisen- 
berg model for both the square-lattice type unfrustrated 
interchain couplings and the triangular-lattice type frus- 
trated interchain couplings. The sublattice magnetiza- 
tion as a function of J2 is then fitted to a power law of 
the form 



M(j2) = (J2 - jry 



(12) 



VI. CONCLUSIONS 

In this paper, we have studied the spin-half and spin- 
one Heisenberg models in spin-chains that arc coupled in 
an anisotropic triangular-lattice geometry, where a spin 
in one chain is coupled to two neighboring spins in an- 
other chain. This geometry is particularly effective in 
preventing the development of spin correlations between 
chains, and in altering spin correlations within a chain. 
For the spin-half case, both Colinear Antifcrromagnctic 
phase and Non-colinear Antifcrromagnctic phase are pos- 
sible in this geometry. Though, in our calculations, the 
Non-colinear phase appears to have a lower energy. Fur- 
ther unbiassed ways of studying this competition numer- 
ically would be useful. Short distance spin correlations 
can be used to distinguish between the different phases. 
Though, potential biassing due to boundary conditions 
should be taken into account. For the spin-one model 
also, we find that this geometry significantly enhances the 
stability of one-dimensional Haldane gap phase, and pre- 
vents the development of long range magnetic order. This 
geometry is quite different from materials like CsNiCls, 
where frustration does not lead to competing correla- 
tions along any given chain and thus the Haldane gap 
phase is quickly destroyed. It would be interesting if ma- 
terials similar to the spin-half materials CS2CUCI4 and 
Cs2CuBr4 are found with spin-one. The study of such 
materials would shed further light on the role of inter- 
chain frustration and the extension of one-dimensional 
physics to higher dimensional systems. 

We would like to thank Ross Mckenzie for discus- 
sions. This work is supported in part by the National 
Science Foundation under grant numbers DMR-0240918 
and PHY05-51164. 



R.R.P Singh: |singh@ raman. physics . ucdavis . edu] ^ S. Sachdev, Low Dimensional Quantum Field Theories for 

A. Auerbach, Interacting Electrons and Quantum Mag- Condensed Matter Physicists (World Scientific, Singapore, 

netism (Springer- Verlag, New York, 1994). 1995), edited by Y. Lu, S. Lundqvist, and G. Morandi. 



7 



^ S. Chakravarty, B. I. Halperin, and D. R. Nelson, Phys. 
Rev. Lett. 60, 1057 (1988). 

* W. Zheng, R. H. McKenzie, and R. R. R Singh, Phys. Rev. 
B 59, 14367 (1999). 

^ B. Bernu, P. Lecheminant, C. Lhuilher, and L. Pierre, 

Phys. Rev. B 50, 10048 (1994). 
® R. R. P. Singh and D. A. Huse, Phys. Rev. Lett. 68, 1766 

(1992). 

Q. F. Zhong and S. SoreUa, Europhys. Lett. 21, 629 (1993). 

* L. Capriotti, A. E. Trumper, and S. SoreUa, Phys. Rev. 
Lett. 82, 3899 (1999). 

^ R. Coldea, D. A. Tennant, R. A. Cowley, D. F. McMorrow, 

B. Dorner, and Z. Tylczynski, Phys. Rev. Lett. 79, 151 
(1997). 

^° R. Coldea, D. A. Tennant, K. Habicht, P. Smeibidl, 

C. Wolters, and Z. Tylczynski, Phys. Rev. Lett. 88, 137203 
(2002). 

R. Coldea, D. A. Tennant, and Z. Tylczynski, Phys. Rev. 
B 68, 134424 (2003). 

R. Coldea, D. A. Tennant, R. A. Cowley, D. F. McMorrow, 

B. Dorner, and Z. Tylczynski, J. Phys.: Cond. Matt. 8, 
7473 (1996). 

M. Y. Veillette, A. J. A. James, and F. H. L. Essler, Phys. 
Rev. B 72, 134429 (2005). 

D. Dalidovich, R. Sknepnek, A. J. Berlinsky, J. Zhang, and 

C. Kallin, Phys. Rev. B 73, 184403 (2006). 

J. O. Fjaerestad, W. Zheng, R. R. P. Singh, R. H. McKen- 
zie, and R. Coldea, Phys. Rev. B 75, 174447 (2007). 

" W. Zheng, J. O. Fjaerestad, R. R. P. Singh, R. H. McKen- 
zie and R. Coldea, Phys. Rev. B 74, 224420 (2006); W. 
Zheng, J. O. Fjaerestad, R. R. P. Singh, R. H. McKenzie 
and R. Coldea, Phys. Rev. Lett. 96, 057201, (2001). 

" W. Zheng, R. R. P. Singh, R. H. McKenzie, and R. Coldea, 
Phys. Rev. B 71, 134422 (2005). 

M. Kohno, O. A. Starykh, and L. Balents, Nature Physics 
3, 790 (2007). 

O. A. Starykh and L. Balents, Phys. Rev. Lett. 98, 077205 



(2007). 

H. J. Schulz, T. A. L. Ziman, and D. Poilblanc, J. Phys. I 
France 6, 675 (1996). 

21 J. Oitmaa and W. Zheng, Phys. Rev. B 54, 3022 (1996). 

22 J. Alicea, O. L Motrunich, and M. P. A. Fisher, Phys. Rev. 
Lett. 95, 247203 (2005). 

2^ J. Alicea, O. I. Motrunich, and M. P. A. Fisher, Phys. Rev. 

B 73, 174430 (2006). 
2'' S. Yunoki and S. SoreUa, Phys. Rev. B 74, 014408 (2006). 
2^ C. H. Chung, J. B. Marston, and R. H. McKenzie, Journal 

of Physics: Condensed Matter 13, 5159 (2001). 
2'' J. Merino, R. H. McKenzie, J. B. Marston, and C. H. 

Chung, Journal of Physics: Condensed Matter 11, 2965 

(1999). 

2^^ M. Q. Weng, D. N. Sheng, Z. Y. Weng, and R. J. BursiU, 

Phys. Rev. B 74, 012407 (2006). 
2^ S. Moukouri, Phys. Rev. B 70, 014403 (2004). 

23 M. Arlego and W. Brenig, Phys. Rev. B 75, 024409 (2007). 
^° M. P. Gelfand and R. R. Singh, Adv. Phys. 49, 93 (2000). 

31 H. Bethe, Z. Phys. 71, 205 (1931). 

32 R. J. Baxter, J. Stat. Phys. 9, 145 (1973). 

33 J. Oitmaa, C. Hamer, and W. Zheng, Series Expansion 
Methods for Strongly Interacting Lattice Models (Cam- 
bridge University Press, New York, 2006). 

3'' M. Bocquet, F. H. L. Essler, A. M. Tsvelik, and A. O. 

Gogolin, Phys. Rev. B 64, 094425 (2001). 
35 L Affleck, J. Phys.: Cond. Matt. 1, 3047 (1989). 
3'5 R. R. P. Singh and M. P. Gelfand, Phys. Rev. Lett. 61, 

2133 (1988). 

3^ S. K. Ma, Modem Theory of Critical Phenomena (Ben- 
jamin Cummings, London, 1976). 

3* L Affleck, Phys. Rev. Lett. 62, 474 (1989). 

33 W. J. L. Buyers, R. M. Morra, R. L. Armstrong, M. J. 
Hogan, P. Gerlach, and K. Hirakawa, Phys. Rev. Lett. 56, 
371 (1986). 



